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SECTION  1 


INTRODUCTION 


There  are  many  problems  of  practical  interest  which  require  a 
mathematical  model  for  a  two-phase  turbulent  flow  of  suspended 
particles  in  a  viscous  fluid.  Most  one-  and  two-equation 
models  for  two-phase  flows  (e.g.,  Refs.  1  and  2)  are  based  on 
ad  hoc  modifications  of  the  single-phase  turbulence-kinet ic - 
energy  and  length-scale  equations  and  fail  to  adequately 
predict  the  physical  behavior  of  two-phase  flows.  Recently, 
Elghobashi  et  al .  (Ref.  3)  presented  a  rigorous  derivation  of 
the  turbulence-energy  and  dissipation-rate  equations  from  the 
momentum  equations  for  an  incompressible  dispersed  two-phase 
flow  and  successfully  predicted  the  main  features  of  a  round 
gaseous  jet  laden  with  uniform-size  solid  particles  (Refs.  4 
and  5).  In  situations  such  as  the  dust  transport  by  turbulence 
in  nuclear  burst  flow  fields  (Ref.  6)  and  supersonic  nozzle 
flows  (Ref.  7)  the  compressibility  of  the  fluid  must  be  taken 
into  account.  The  objective  of  the  present  work  is  to  extend 
the  two-equation  turbulence  model  in  Ref.  j  to  become  appli¬ 
cable  to  a  compressible  dispersed  two-phase  flow  in  which  the 
fluid  has  variable  density  and  the  particles  have  constant 
density.  The  extended  model  can  also  be  used  to  study  such 
problems  as  compressible  adiabatic  mixing  and  low-speed  iso¬ 
thermal  mixing  of  two  dissimilar  two-phase  flows  (Ref.  8). 
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SECTION  2 


EQUATIONS  OP  MOTION 

We  begin  the  formulation  of  the  problem  by  stating  the 
assumptions  involved  in  deriving  the  equations.  These  are: 

1.  Both  phases  behave  macroscopically  as  a  continuum, 
but  only  the  carrier  fluid  behaves  microscopically 
as  a  continuum  with  variable  density.  This  means 
that  the  volume-averaged  equations  are  based  on  a 
control  volume  larger  than  the  particle  spacing  but 
much  smaller  than  the  characteristic  volume  of  the 
flow  system.  Mutual  exclusion  of  the  phases  is 
also  ensured. 

2.  The  dispersed  phase  consists  of  rigid  particles 
spherical  in  shape,  uniform  in  size  and  constant  in 
density.  The  uniformity  of  size  reduces  the 
magnitude  of  bookkeeping  at  this  stage  of  the  work, 
and  thus  concentrates  the  effort  on  understanding 
the  mechanisms  of  interactions  between  the  two 
phases.  Extension  to  nonuniform  size  distribution 
is  a  straightforward  matter  (Ref  9). 

3.  The  volume  fraction  of  the  dispersed  phase  is 
such  that  no  collisions  occur  between  the 
particles.  This  assumption  renders  the  equations 
valid  only  for  dilute  suspensions. 

4.  Neither  the  suspended  matter  nor  the  carrier  fluid 
undergoes  any  phase  changes.  Although  this 
assumption  rules  out  some  situations  of  practical 
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5. 


interest,  it  is  necessary  to  investigate 
complexities  in  a  stepwise  manner. 


Additional  assumptions  on  modeling  of  some  of  the 
turbulent  correlations  are  stated  in  Sections  4, 
and  6.  The  sparseness  of  experimental  data  for 
variable  density  flows  allows  us  to  assume  that 
the  forms  and  values  of  the  coefficients  of  the 
turbulent  correlations  of  constant-density  flows 
apply  to  variable-density  flows  as  well. 


The  instantaneous,  volume-averaged  momentum  equations,  in 
Cartesian  tensor  notations,  of  the  carrier  fluid  are  thus 
(Refs.  3  and  10) 


( Qlui  >  ,  t  +  (QlUjUi>,j  =  -<l-K+2>P,i  -  KF+^Ui-Vi) 

+  [Mi+itUi.J+Oj.i)],  j 

-  i(Hi+1Uj,J|),i  +  QlSfi  +  Fli  • 


The  corresponding  equations  for  the  particle  phase  are 
<Q2vi>,t  +  (Q2vjvi>,j  =  -  +  2p,i  +  F^tUi-Vi) 

+  [^2+2  <  vi  ,  j+vj  ,  i  ^  ]  ,  j  “  §(H2*2vjJ,  Jl5  ,  i 


'W-.W-WsWs: 


1 

V 

V 


$ 


*■;  -rf-;  v-.'v.;  j-\  -f  1  VJ  V";  IP 


The  contir.uitv  ecuaticns  are 


Q  1  t  f  ( Q:U;  )  (  i  =  0 


for  the  fluid  and 


-  -  (  Q  o  V  j  l  .■  =  3 


for  the  Darticle  chase.  The  alobai  continuitv  ecuati. 


?  +  D  =  j_  =  £l  + 

1  2  "  "  °1  ~2 


(5) 


In  Eqs .  (l)  -  (5)  and  throughout  the  report,  the  subscripts 

1  and  2  denote  the  fluid  and  particle  phase  respectively. 
Partial  derivatives  are  represented  by  a  subscript  consisting 
of  a  comma  and  an  index  [e.g.  ,  (  )  «.  =  3i  )  5t:  L'j  ,  s  IT .  3::. 

U  <  k;J  -  3  i  3xk3xjJ  where  K;  fi  =  1,2.3)  are  the  rectangular 
spatial  coordinates],  are  the  velocity  components  of  -he 
fluid,  7j  are  the  velocity  components  of  the  particle  phase .  o 
and  u  are  the  material  density  and  viscosity.  J.  =«■•>;  os  the 
apparent  density,  <*>  is  the  volume  fraction.  P  is  the  press'.:-'  . 
g<  is  the  component  of  gravitational  acceleration  in  the  : 
direction.  F;  is  the  component  of  body  forces  other  than  •"hat 
due  to  gravity,  and  F  is  the  interphase  friction  coeffi -oent 
=  1  ;  Mi  d-  for  Stokes'  flow  around  a  gar~i:.e  o:  ;!  i  amo  *  -  1 

K  is  toe  local  effectiveness  of  momentum  transfer  fr:.r.  -  e 
particle  phase  to  the  fluid  and  is  discussed  in  dt'-ai  i 
Per.  II.  It  suffices  here  to  state  to.at  K  equals  uni~y  for 
accelerated  particle  phase  and  assumes  lower  values  when  - 
phase  is  decelerated:  the  minimum  value  of  K  being  cere.  Ic. 
teneral .  d  deoends  on  the  local  orooertias  of  the  : 1  old  and 


! 


3 


3 

m  * 


a 


a 


k 


'vJ 


turbulence,  the  slip  velocity,  and  the  particulate  size  and 
concentration . 

The  mean  flow  equations  are  now  obtained  from  the 
instantaneous  ones,  for  variable  ,  constant  p2  and  n-j.  and 
zero  ^2  (consistent  with  the  dilute  suspension  approximation 
stated  in  Ref.  10,  p.  256),  by  performing  the  conventional 
Reynolds  time-averaging  of  Eqs .  (1)  to  (5).  (The  density- 

weighted  averaging  of  Favre,  Ref.  11,  does  not  render 
significant  simplification.  We  use  Reynolds  averaging  mainly 
because  most  experimental  data  refer  to  time  averaging 
correlations).  The  mean  momentum  equations  of  the  fluid  are 

(QlUi  +  q^7),t  +  (Qluj°i>  ,  j  =  ~  (l-K+2)P,i  +  K*2p(i 

-  KF^<fr2(Ui-Vi)  +  j] 

+  Ml[*1(01>:j+Uj>1)  +  *i(uif  j+uj#1)]f  j 

-  §  Hi<+i ,i 

-  (QiUiUj  +  Oi  qxUj  +  Uj  q1ui  + 

+  Qigi  +  F1A  .  (6) 


5 


The  mean  momentum  equations  of  the  particle  phase  are 


‘WVi’.t  +  (Q2VjV,j  '  ->»2P,i+«2P,i1 


+  F(t2(Ui-Vi)  +  ♦2(ui-vi)] 


-  (Q2viVj  +  Viq2vj  +  q2vivj 


+  ®2®i  +  p2i  • 

The  mean  continuity  equation  of  the  fluid  is 

Qi.t  +  <Qiui  +  qiuiT,i  =  o 

The  mean  continuity  equation  of  the  particle  phase  is 


Q  2  ,  t  +  (  Q  2^  i  +  *5  2V  i  )  ,  i  =  0 
which  can  be  written  as 


*2  ,  t  +  (  +2vi  +  *2vi>  ,  i  =  0 

since  p 2  is  constant.  The  mean  global  continuity  equati 

+i  +  +2  =  1 

which,  when  substrated  from  Eq .  (5),  gives 

$  1  +  $2  =  0 


t  j 

(7) 

(8) 

(9) 

(  10  ) 

on  is 

(  11  ) 

(  12  ) 


In  Eqs .  ( 6 ) - ( 12)  capital  letters  (except  K  and  F)  denote 

time-mean  quantities,  lower-case  letters  (except  pj  and  g* ) 
designate  fluctuating  components  and  overbars  indicate 


Reynolds-averaged  correlations.  For  constant  quantities 
p 2  and  9i)  the  mean  and  instantaneous  values  are  equal,  whereas 
for  all  other  variables  the  instantaneous  value  consists  of  a 


mean  component  and  a  fluctuating  component  (e.g.,  Q  +  q,  U  +  u. 


Pi  +  Pi 


etc .  ) 


SECTION  3 


TURBULENCE-KINETIC-ENERGY  AND  DISSIPATION-RATE  EQUATIONS 

The  first  step  in  the  derivation  of  the  equations  of  the 
fluid's  turbulence  kinetic  energy  (k  s  uiu^/2)  and  its 
dissipation  rate  {«  s  ^  u^^  where  uj  =  is  to 

obtain  a  transport  equation  for  by  subtracting  Eq.  (6) 
from  Eq.  (1).  The  k  equation  is  produced  by  multiplying  the 
Uj  equation  throughout  by  u^  and  then  time-averaging.  The  « 
equation  is  obtained  by  differentiating  the  u^  equation  with 
respect  to  xk,  multiplying  throughout  by  2i/1ui  lc  and  finally 
time-averaging . 


The  resulting  k  and  €  equations  are  given  in  Appendix  A. 
The  closure  of  these  equations  is  discussed  in  Sections  5 
and  6,  respectively. 


SECTION  4 


CLOSURE  OF  THE  MOMENTUM  EQUATIONS 

The  turbulent  correlations  appearing  in  Eqs .  (6)  and  (7)  are 

of  five  types: 

1.  Correlation  of  velocity  fluctuations  with  those  of 
the  volume  fraction  or  apparent  density,  e.g., 

or 

2.  The  pressure  interaction  correlation  <^P 

3.  Multiple  correlations  among  various  components  of 
velocity  fluctuations  with  those  of  the  apparent 
density,  e.g.,  q^jUj  ; 

4.  Correlations  of  strain  rate  fluctuations  with  those 
of  the  volume  fraction,  e.g.,  ^1uiij  ; 

5 .  Multiple  correlations  of  various  components  of 
velocity  fluctuations,  e.g.,  u^uj  . 

The  first  four  types  occur  only  due  to  the  presence  of  the 
second  phase;  their  modeling  is  discussed  below. 

According  to  Ref.  3  ( Eq .  (10))  and  Ref.  12,  we  model  the 
turbulent  flux  by  a  gradient  transport  term  and  a 

convective  transport  term  such  that 


where  vt  is  the  kinematic  eddy  viscosity  (=  c^  '<z  '  e  with  c^  = 
0.09,  Ref.  13,  under  the  provisional  assumption  that  it  has  the 
same  value  as  in  the  case  of  constant  density)  and  a ^  is  the 

turbulent  Schmidt  number  of  Similarly  we  model  q^Uj.  by 

q1ui  *  * < ^ Q1 , i  "  2  {vt/ffq),i  (14) 

where  tr^  is  the  turbulent  Schmidt  number  of  q.  As  long  as 
experimental  data  for  are  not  available  we  make  a 
provisional  assumption  that  «r  is  equal  to  «r^  (=  1.0  for  the 
sample  calculation  in  Ref.  3). 

According  to  Ref.  3  ( Eqs .  (11),  (12),  (14),  (15)  and  (16))  and 

Refs.  14,  15  and  16  we  model  <£2p  <  by 
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where  ea  and  en  are  unit  vectors  and  the  approximate  values  of 
the  coefficients  are 
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According  to  Ref.  16  ( Eq .  (6.45))  we  model  q^^Uj  by 

qiuiuj  ■  "c$5  (k/e  )  [  viiUjj  (ujqi)  #  jj  +  UjUj,  (u^)^]  ,  (17) 

where  the  proportionality  constant  c^5  is  approximately 
equal  to  0.1. 

The  strain-rate  volume-fraction  correlations  of  the  type 
^lui,j  only  appear  multiplied  by  the  molecular  viscosity  of 
the  fluid  and  therefore  will  be  neglected  due  to  its  relatively 
small  magnitude. 

The  last  correlation  to  be  modeled  in  the  momentum  equations 
is  that  of  the  form  u^Uj.  Again,  to  be  consistent  with  the 
present  level  of  closure,  this  quantity  will  be  calculated 
from  (Eq.  (18)  in  Ref.  3) 

uiu  j  =  "  vt ( ui ,  j+uj  ,  i  )  +  3  sijumum  ~  3  vt  *ijuj,j  •  (18) 

This  completes  the  modeling  of  the  momentum  equations. 


SECTION  5 


CLOSURE  OF  THE  TURBULENCE  KINETIC-ENERGY  EQUATION 

The  exact  equation  of  the  turbulence  energy  k  for  the  carrier 
fluid  appears  in  Appendix  A  and  consists  of  34  terms.  They 
are  classified  into  groups  enclosed  by  square  or  curly 
brackets;  each  group  is  labeled  according  to  its  particular 
contribution  to  the  conservation  of  k. 

The  various  correlations  in  these  groups  range  from  second  to 
fourth  order.  We  decide  at  the  outset  on  neglecting  all 
fourth-order  correlations  such  as  <31^ (uj.u  j )  (  j  and  UjUjUjq1;j 
Also,  the  contribution  to  the  diffusion  of  turbulence  energy 
due  to  the  pressure  interaction  (u^p)  ^  will  be  neglected  as 
it  is  of  relatively  small  magnitude  (Ref.  17).  Now  the 
remaining  terms  will  be  modeled. 

The  five  transient  terms  will  be  collectively  approximated  by 
(  Ql lc )  #  t  *  The  convection  terms  require  no  approximation.  The 
production  group  contains  the  correlations  q^^uj  and  u^u-f 
which  have  been  evaluated  earlier  by  Eqs .  (17)  and  (18).  The 

pressure  velocity-divergence  correlation  pu^ ( ^  in  the 
turbulent-diffusion  group  cannot  be  neglected  here  since  u^  i 
does  not  vanish  in  two-phase  flows.  pu^  ^  is  evaluated 
following  the  approach  of  Ref.  18,  thus  (Eqs.  (19)  and  (20)  o 
Ref.  3)) 
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where 


pii  =  -2(uiu:<  ui,k: 


p  *  1  pii 


*ii  =  “2(uiuk  Uk,i> 


c:  =  1.5  ,  c  2  =  0.4 


The  last  two  terms  In  the  turbulent-diffusion  group  can  be 
modeled  as  ( Eqs .  (21)  to  (25)  and  Eq.  (10)  of  Ref.  3) 


uiui(Qiuj),j  *  Qiuiujui,j  =  Qih/ffk)k,j 


+  2  k<lVk,.j  ,j 


where  the  turbulent  Schmidt  number  of  k  is  taken  as  or  =  i.o 
(Ref.  13,  under  the  provisional  assumption  that  it  has  the 
same  value  as  in  the  case  of  constant  density) . 


There  are  eight  terms  in  the  extra  production  and  transfer 
group,  the  last  two  of  which  are  neglected  for  being  of  fourth 
order.  The  remaining  six  terms  are  modeled  next. 


The  second  term,  -Q.U.  u-u.  is  modeled  following  the 

1  i  1  J  /  J 

proposal  of  Ref.  19  as  Q-U.1/..U..  ,  ..  The  correlation  of  the 
_  lit  j,ij 

form  u.(q.u.)  .  which  appears  in  the  third  and  fourth  terms  is 

1  +  J  f  J 

expanded  as 


where  the  first  term  on  the  right  is  evaluated  using  Sq.  (17) 
and  the  second  term  is  modeled  as 
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where  the  provisional  assumption  on  has  already  been 
mentioned  in  Section  4.  The  correlations  (in  the  first 

term)  and  q^u^  (in  the  fifth  term)  have  been  discussed  earlie 


(Eqs.  (13)  and  (14! 


The  sixth  term  is  approximated  by 


where  q1ui  j  is  neglected  for  being  relatively  smaller  than 

(’4IqT>,j  : 

The  extra  dissipation  group  contains  three  terms  which  exist 
only  due  to  the  slip  between  the  two  phases..  According  to 
Ref.  3  (Eq.  (39))  and  Refs.  20  and  21  we  model  the  correlation 
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in  which 


12u1/p1d2 

(31  ) 

3p1/(2p2  +  p,) 

;  32 

and  V  is  Taylor's  microscale.  The  last  term  in  the  extra 
dissipation  group  contains  the  triple  correlations  ^u.u,  and 

^2uivi  which  can  be  modeled  by  Sqs .  (17)  and  (42)  of  Ref.  3 

as 
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and 
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where  the  double  correlations  on  the  right  sides  have  been 
modeled  earlier  ( Eqs .  (13)  and  (13)). 


The  term  >i,  +  ,  u<(u4  .,+u.  .  )  ..  constitutes  the  dissipation  of 

due  to  viscous  action,  and  if  ♦.  is  set  to  unity  it  : 


is  reduce 


to  the  single-phase  dissipation  terms.  This  term  is  modeled  as 
(Eq.  (43)  in  Ref.  3) 


*1+1  Ui(ui<;J+ujfi)  3  -#»!*!€  (35) 

The  six  terms  in  the  viscous  diffusion  and  dissipation  group 
will  be  neglected  due  to  their  relatively  small  magnitudes  as 
compared  to  the  terms  in  the  turbulent  diffusion  group. 

The  first  term  in  the  field  forces  effects  group  is  modeled  by 
Eq.  (14).  Similarly  the  second  term  in  the  group  can  be 
modeled  as 


uifli  -  )Fllf  i  -  2*  n  (vt/flrf  )  ,i 


(36) 


where  we  make  a  provisional  assumption  that  9^  equal  to  or^  . 


SECTION  6 


CLOSURE  OF  THE  TURBULENCE-ENERGY  DISSIPATION-RATE  EQUATION 

The  exact  equation  of  the  dissipation  rate  of  turbulence 
energy  e  for  the  carrier  fluid  appears  in  Appendix  A  and 
consists  of  52  terms.  They  are  classified  into  groups  similar 
to  those  of  the  k  equation. 

Again  we  neglect  all  fourth-order  correlations  as  mentioned 
in  the  previous  section. 

All  the  terms  in  the  first  group,  labeled  Transient,  are 
approximated  collectively  by  (Q^e)  t  . 

The  convection  group  consists  of  eight  terms  of  which  only 
the  first  and  the  second  are  of  higher  magnitude  than  the 
other  six  at  large  Reynolds  number.  This  is  based  on  an  order 
of  magnitude  analysis  (Ref.  22)  which  shows  that  the  first  and 
second  terms  are  greater  than  the  others  by  at  least  a  factor 
of  (J/X),  which  is  of  order  (Rjj)1/^  Here  J3  is  the  length 
scale  of  the  energy  containing  eddies,  X  is  a  Taylor’s 
microscale,  and  R^  is  the  Reynolds  number  based  on  J? . 

The  third  term  in  the  production  group  is  decomposed  as 


-2w1Q1(uiUjfku1#k) (j  =  -2v1Q1(uif jUj(kui(k  *  uiUj|kuik 


+  uiuj ,kui,kj>  '37' 

The  third  term  on  the  right  side  of  Eq .  (37)  and  the  second 

term  in  the  production  group  differ  only  in  their  signs  and 
thus  cancel  each  other.  The  first  term  on  the  right  side  of 


Eq.  (37),  which  represents  the  production  of  e  by  self¬ 
stretching  of  the  vortex  tubes,  is  the  dominant  one  at  large 
Reynolds  number.  It  is  larger  than  the  second  term  by  a  factor 
of  Rj  and  larger  than  the  first  terra  in  the  production  group 
by  a  factor  of  R . 

We,  therefore,  retain  only  ~2v^Q1  ui ,  ju  j  ,  kui  ,  k  as  t^le  ma-in 
generation  of  «.  This  term  and  the  extra  production  terms  are 
modeled  collectively  as 


-2v1Q1  ui(jUj  kuik  +  extra  production  terms  =  c€lGk«/k  ,  (38) 

where  Gk  is  the  total  production  of  k  discussed  in  Section  5, 
and  c€l  is  a  constant  of  value  1.43  (Ref.  13,  under  the  provi¬ 
sional  assumption  that  it  has  the  same  value  as  in  the  case  of 
constant  density) .  "Total"  here  means  the  production  terms  which 
are  common  to  the  single-phase  and  two-phase  k  equations  in 
addition  to  the  extra  production  and  transfer  terms. 

The  turbulent  diffusion  group  contains  six  terms.  At  high 
Reynolds  number  only  the  last  two  terms  will  be  retained;  they 
are  larger  by  at  least  a  factor  of  Rjj/2  than  the  other  terms. 
These  two  terms  will  be  modeled  collectively  as 

-2w1(Q1ujuirkuif +  Qx , j  ujuifkulfk) 

=  [Qj  (vt/«r€  )€  (  j  +  \  Qle  (vt/9€  )  ,  j]  ,  j  (39) 


All  the  terms  in  the  extra  production  group  except  the  mean 
pressure  gradient  term  are  smaller  than  the  main  production 
term,  modeled  in  Eq .  (38),  by  at  least  a  factor  of  R^1  2  and 


thus  can  be  neglected.  The  mean  pressure  gradient  term  is 
included  in  the  production  term  { Eq .  (38)). 

The  first  term  in  the  viscous  diffusion  and  dissipation  group 
represents  the  main  dissipation  of  e;  it  reduces  to  the 
single-phase  form  when  <l»1  equals  unity.  This  term  is  larger 
than  the  other  terms  in  the  group  by  a  factor  of  R^/2  and  thus 
it  is  the  only  one  retained.  Now  the  total  dissipation  of  e 
includes  this  term  in  addition  to  the  extra  dissipation  terms. 
They  are  modeled  collectively  as  Q1(€/k)(c€2«  +  c€3€e)  where 
«e  represents  the  extra  dissipation  terms  appearing  in  the  k 
equation,  cfi2  is  a  constant  of  a  value  about  1.92  (Ref.  13)  and 
c^3  is  a  constant  of  value  1.2  (Ref.  3)  (both  under  the 
provisional  assumption  that  they  have  the  same  values  as  in 
the  case  of  constant  density) . 

Both  terms  in  the  field  forces  effects  group  appear  multi¬ 
plied  by  the  kinematic  molecular  viscosity  of  the  fluid  and 
therefore  will  be  neglected  due  to  their  relatively  small 
magnitudes . 


SECTION  7 


A  SAMPLE  APPLICATION 

As  an  example  of  the  application  of  the  modeled  k  and  e  equa¬ 
tions,  let  us  consider  the  motion  of  the  dusty  air  during  a 
nuclear  explosion  (Ref.  6).  In  order  to  understand  the  essen¬ 
tial  features  of  the  complicated  phenomena,  some  highly  ideal¬ 
ized  and  well-controlled  laboratory  tests  have  been  performed 
or  planned  such  as  high  speed  wind  above  a  sand  bed  in  a  wind 
tunnel  (Ref.  23)  and  shock  wave  sweeping  a  sand  bed  in  a  shock 
tube  (Ref.  24).  To  interpret  and  correlate  the  results  from 
such  tests,  we  may  use  the  two-dimensional  version  of  the 
present  model  in  Cartesian  coordinates  x  and  y.  The  final 
form  of  the  modeled  k  and  e  equations  is  given  in  Appendix  B 
(under  the  assumptions  that  the  diffusional  fluxes  in  y  direc¬ 
tion  are  much  larger  than  those  in  x  direction,  that  the 
effective  coefficient  K  *  1  and  that  gy  =  g  and  gx  =  Flx  =  Fly 
=  0)  ,  The  remaining  modeled  equations  contain  two  mean  con¬ 
tinuity  equations  (from  Eqs .  (8)  and  (9)),  one  mean  global 
continuity  equation  (Eq.  (11)),  four  mean  momentum  equations 
(from  Eqs.  (6)  and  (7)  and  Section  4)  and  one  mean  equation 
of  state  such  as  the  perfect  gas  equation 

P  =  RPiTi  (40) 

where  R  is  the  specific  gas  constant  and  Tj_  is  the  absolute 
temperature  of  the  air.  Thus,  for  isothermal  problems  (T^  = 
constant)  we  have  ten  equations  for  ten  unknowns  Q^,  q2,  . 

Ux,  Uy,  Vx,  Vy ,  P,  k  and  e  (from  which  =  Q1/p1  and  +2  = 
Q2/p2  can  t,e  readily  obtained)  .  With  proper  initial  and 
boundary  conditions,  these  equations  can  be  solved  numerically 
by  a  marching  finite-difference  procedure  described  in  Ref.  5. 


The  results  of  such  calculation  will  be  presented  in  a  forth¬ 
coming  report. 

If  is  variable  and  the  energy  exchange  between  the  air  and 
dust  can  be  neglected  we  need  to  include  a  mean  energy  equa¬ 
tion  for  the  air  in  the  numerical  solution  procedure. 

If  the  energy  exchange  between  the  air  and  dust  (at  tempera¬ 
ture  T2)  is  not  negligible,  we  need  to  include  the  mean  energy 
equations  for  both  air  and  dust  in  the  numerical  solution. 


SECTION  8 


CONCLUSION 

The  k  -  «  turbulence  model  for  an  incompressible  dilute 
suspension  of  Ref.  3  has  been  extended  to  a  compressible 
dispersed  two-phase  flow  by  introducing  the  apparent  densi¬ 
ties  and  Q2  and  the  material  density  as  new  variables. 
The  fluid  has  variable  density  and  the  particles  have  con¬ 
stant  density.  This  allows  the  application  of  the  model  to  a 
wider  class  of  practical  problems. 

As  in  Ref.  3,  the  k  and  e  equations  are  first  rigorously 
derived  from  the  two-phase  momentum  equations  and  then  their 
closure  is  provided.  This  is  in  contrast  to  the  usual 
approach  based  on  ad  hoc  modifications  of  the  single-phase 
turbulence-kinetic-energy  and  length-scale  equations. 

The  proposed  closure  of  the  equations  accounts  for  the  inter¬ 
action  between  the  two  phases  and  its  influence  on  the  turbu¬ 
lence  structure.  Sparseness  of  experimental  data  for  variable 
density  flows  necessitates  some  provisional  assumptions  that 
forms  and  values  of  the  coefficients  in  the  turbulent  corre¬ 
lations  of  constant-density  flows  apply  to  those  of  variable- 
density  flows.  Such  assumptions  indicate  areas  of  needed 
experimental  investigations  which,  when  completed,  can  in  turn 
modify  the  present  work  and  enhance  its  validity. 
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APPENDIX  A 


EXACT  EQUATIONS  OF  KINETIC  ENERGY  OF  TURBULENCE 
AND  DISSIPATION  RATE  OF  THAT  ENERGY 


The  exact  equation  of  the  turbulence  kinetic  energy  (k  =  UjUj/2) 
of  the  carrier  fluid  is 

[( Qi  UjUj/2)  (t  +  Q1)t  UiUj/2  +  ui,t  ui<3l  +  ui  uicJl,t 

Transient 

+  ui(<Jl^i),t  ]  +  [(QiUjk)  +  ktQiUj)  (j]  = 

Convection 

-  [( Q1Ui  )  (  j  uiU:j  +  U^jqjUiUj  +  Uj/jqiUiu7] 

Production 


(l-K+2)  [(uAp) , i-puif i]+K^iUiPfi  +  uiui(Q1uj)  (  j 

Turbulent  Diffusion 


+  QiUiUjUi^j  -  [*P,i*lui  +  Qlui  uiuj , j  +  ui 


+  Uj  UitqjUi)^  +  (UiUjJ.jqiUi  +  U^jU^^ 

Extra  Production  and  Transfer 


+  q1ui(uiu;j)  #  j  +  uiUiUjqlfj]  +  KF  [+2 


+  (Vi-Ui)^2ui  +  ♦^itVi-Ui)]  +  ii1[*1  ui(uij+uj(i)fj] 
Extra  Dissipation  Dissipation 
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rREVIOUSP*GE 

IS  BLANK 


Hi  +1,3  ui<Ui,j 


i.  jTUJ.i 


>*l],jui  "  l[(  +  luJ>,  J»>  ,  1U1  +  ^1ujJ,jJ) 
Viscous  Diffusion  and  Dissipation 


+  ^1uj»,j»)  ,  iui J  |+  1 9 i  ui9 1  +  Mfli 


Field  Forces  Effects 


The  exact  equation  of  the  dissipation  rate  of  k 
5  V1  ui,kui,k>  is 

^[(QiUi)  ,ktuifk  +  (qi«i)  ,ktui,k  +  (9iUi)  ,ictui,k] 

Transient 

+  |  ( Q  iU  j «  )  ,  j  +  *(QlUj),j  +  2w1[(Q1Ui) 

+  (Ql°i)  ,kuj, jui,k  +  (Qlui)  , jkui,kuj  +  ( Qlu j  >  ,kui,kui, j 

Convection 


+  (Q 


luj>,jk  uiui,k  +  Qlui  u i , ku j , j k] j  =  ~2vl [Q1  ui,kui,kui 


-  ^1  uiui,jkuj,k  +  Ql(uiuj,kui,k) , j] 


Production 


+  2v1KF| [+2(vi-ui)j (kuiik  +  [^2<vi-ui>] ,k^i, 


Extra  Dissipation 


+  U2(Vi-Ui)]<kui(k  [-  2v1[Q1(Jk  uiUjuifk 


+  Qi,k(uiuj,  j+ujui,  j>ui,k  +  Qi.j  uiuj,kui,k  +  P,ki’Ji,k 

Turbulent  Diffusion 


+  Qlujui,kui, jk  +  Ql, jujui,kui,k]  +  2w1j(K*2?, 


,  kai  ,  .< 


~  KP,ki  ^  lui  ,  k  "  p-i  ui,k<K4>l),k  ~  <  K<*  lP  ,  i  >  ,  ku  i  ,  k 


"  ui,k[<3lui>  ,  jaj  ,k+^lui)  ,kuj  ,  j]  -  ui  3lui,kuj.;l 


u  jui  ,  k(  *3 1^  i  )  ,  jk  "  uj  £5iui,kui,jk  "  ui  ( 3  lu  j  *  ,  jkui  ,  k 

Extra  - 


^  ( <5 1 U j  )  "  ju^  (  kui  ,  k  +  ^luj  ^  ,  kui  ,  jui  ,  k] 

-  Production 

[uiuj^l ,  jk^l  ,k<[Jiuj  ,  j+ujui  ,  j  )]ui  ,k 
3i  (  ui  ,  juj  )  ,  k+uiuj  ,  k^l ,  j  +c*i  (  uj  ,  j  ui  )  ,  k]  ui  ,  k 


-  ujui,kSi, jui(k  -  [ui^j<?i, jk+<3i,k<uiuj,:+ujui, j>]ui,k 

-  5i  (  UiU  J  f  k  )  f  j  -MliUj  ,  kqi  ,  j  +qi  (  u  j  Ui  f  k  )  (  j  +U-j  Ui  (  kqi  ,  j  ]  Ui  ,  k" 


+  2|Dll/lJ  [*l(ui,  j+uj,i>]  ,  jkui,k  +  [^i(Ui(  j+Uj(i)](,kui( 


[^1  (  ui  ,  j +u  j  ,  i  1  ]  ,  jkui  ,  k  “  3  [(  +  luJJ,  jJ  >  ,  ikui  , 
Viscous  Diffusion  and  Dissipation 


i ,  k 


+  uj? ,  A  )  ,  ikui  ,  k  +  (  ^  1UJ? ,  j5  )  ,  ikui  ,  k 


+  2vl  [9iqi,kui,k  +  fli,kui,kF  • 
Field  Forces  Effects 
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APPENDIX  B 

THE  MODELED  FORM  OF  <  AND  e  EQUATIONS 

The  modeled  conservation  equations  of  the  kinetic  energy  and 
the  dissipation  rate  of  that  energy  for  the  carrier  fluid  in 
the  sample  application  are  listed  here. 

(i)  The  k  equation: 


Production  (P) 


Extra  Production  (P  ) 

-i  e 
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\p2/[ 


*\  J0  “2  - /  ry  Y  ly  ^  J  2  ,y 

Extra  Dissipation  (e  ) 


Dissipation  Field  Forces  Effects 


(ii)  The  e  equation: 


Q,  £  . 

+  Q.  U  e  +  U  £ 

1  ,t 

l  x  ,x  y  ,y 

Transient  Convection 


-  C£l  <£/*)(  P+Pe)|  +  [°l(Va£Ky  +  €(vt/ac),: 


Total  Production 


Diffusion 


-  Q^c/k)  (ce2  e  +  ce3  £e) 
Total  Dissipation 


The  notations  used  in  the  partial  derivatives  have  been 
explained  in  Section  2  of  the  text,  thus  (  )  means 

t  ^ 

3(  )/3t,  (  )  means  3(  )/3x  and  (  )  means  3(  )/3y,  where 

' x  i  y 

x  and  y  are  the  distances  along  the  horizontal  and  vertical 
directions,  respectively.  The  values  of  the  constants 
appearing  in  the  two  equations  are : 


c*  1 


=  1.43,  c  =  1.92,  c  =  1.2,  =  0.1,  crv  =  1.0 
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APPENDIX  C 


NOMENCLATURE 

H'  C€l'  °€2'  Ce3'  C$l'  C4>2'  C  <p3'  C  <p4'  C<p5 
constants  in  the  turbulence  model 
particle  diameter 

an  expression  defined  in  Eqs .  (20) 

unit  vectors 

an  expression  defined  in  Eq.  (30) 

fluctuating  component  of  body  forces  other  than 

that  due  to  gravity 

instantaneous  (in  Eqs.  (1)  and  (2))  or  time-mean 
component  of  body  forces  other  than  that  due  to 

gravity 

interface  friction  coefficient 
component  of  gravitational  acceleration 
total  production  of  k 

turbulence  kinetic  energy  of  the  fluid 

local  effectiveness  of  momentum  transfer  from  the 

particle  phase  to  the  fluid 

length  scale 

pressure  fluctuation 

instantaneous  (in  Eqs.  (1)  and  (2))  or  time-mean 
pressure 

an  expression  defined  in  Eqs.  (20) 
production  (in  Appendix  B) 
extra  production  (in  Appendix  3) 
fluctuation  of  the  apparent  density 
instantaneous  (in  Eqs.  (  1 )  —  ( 5 ) )  or  time-mean 

apparent  density 
specific  gas  constant 
Reynolds  number  based  on  Jl 


T  absolute  temperature 

Ui  fluctuating  velocity  component  of  the  fluid 

ui  instantaneous  (in  Eqs .  ( 1 )  —  ( 3 ) )  or  time-mean 

velocity  component  of  the  fluid 

vi  fluctuating  velocity  component  of  the  particle 

phase 

vi  instantaneous  (in  Eqs.  (1)— (4) )  or  time-mean 

velocity  component  of  the  particle  phase 

xi  rectangular  spatial  coordinate 

x  horizontal  coordinate 

y  vertical  coordinate 

Greek  symbols 

*  an  expression  defined  in  Eq.  (31) 

(l  an  expression  defined  in  Eq.  (32) 

Kronecker  symbol 
«  dissipation  rate  of  k 

€e  extra  dissipation  terms  in  the  k  equation 

X  Taylor's  microscale 

viscosity 

v  kinematic  viscosity 

wt  kinematic  eddy  viscosity 

p  material  density 

«</>.  «fq«  <*k»  °-f  ,  *€ 

turbulent  Schmidt  numbers 

$  fluctuation  of  the  volume  fraction 

♦  instantaneous  (in  Eqs.  (1)— (5) )  or  time-mean  volume 
fraction 

w  frequency  of  turbulence 

an  expression  defined  in  Eq.  (27) 

an  expression  defined  in  Eq.  (28) 

Ar  an  expression  defined  in  Eq.  (29) 


I 


I 


Subscripts 

1  fluid  phase 

2  particle  phase 

,t  partial  derivative  with  respect  to  t 

,  j  partial  derivative  with  respect  to  Xj 

Superscript 

-  time-averaged  value 

'  fluctuating  component 
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